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1.  statement  of  the  Problem 


1.1  Background 

The  main  objective  of  this  project  is  to  determine  if  heat  pipe  technology  can  be 
used  to  cool  the  portion  of  a  missile's  control  surfaces  heated  by  hot  exhaust  gas.  The 
goal  is  to  reduce  the  heat  conducted  into  the  body  of  the  missile  and  its  actuation  rod. 
The  small  portion  of  the  airfoil  near  the  missile  body  heated  by  the  exhaust  plume  is  to 
be  cooled  by  rejecting  the  heat  to  the  large  portion  of  the  fin  moving  through  ambient 
air.  It  is  a  ground-to-ground  missile  with  a  maximum  speed  of  about  Mach  0.75,  but 
the  evaporator  will  see  exhaust  gas  (after  shock)  at  Mach  0.89  and  1213  K.  Figure  1  is 
a  simplified  illustration  of  the  physical  situation. 


Figure  1.  Simplified  conceptual  drawing  of  the  missile  configuration. 

(Note:  There  are  also  two  vertical  fins,  one  above  and  one  below  the  missile  body) 


The  ambient  air  temperature  is  assumed  to  be  300  K,  but  since  this  missile  could  be 
deployed  anywhere,  temperatures  as  low  as  240  K  are  possible  during  transport  or 
even  operation.  Thus,  a  frozen  startup  scenario  is  possible  for  even  room  temperature 
working  fluids. 

1.2  Physical  Requirements 

The  heat  pipe  is  to  be  built  inside  a  NACA  0015  airfoil  of  chord  length  0.1143  m 
(see  Figure  2).  The  total  length  of  the  heat  pipe  will  only  be  a  small  portion  of  the 
0.254  m  span  of  the  fin.  The  first  0.05  m  of  the  fin  will  be  heated  by  the  exhaust  plume. 
The  heat  pipe  shell  has  to  withstand  a  pressure  difference  of  approximately  1.7 
atmospheres  during  flight.  The  heat  pipe  will  be  subjected  to  a  launch  acceleration 
load  of  2g  along  the  launch  flight  path  until  a  maximum  speed  of  Mach  0.75  is  reached 
This  acceleration  is  expected  to  be  short  lived  (~20  seconds).  The  bottom  fin  on  the 
missile  will  be  placed  in  an  evaporator  up  configuration  requiring  the  working  fluid  to 
wick  up  against  gravity.  Little  maneuvering  is  anticipated  until  the  final  few  seconds  of 
flight.  There  must  be  an  axial  hole,  0.0159  m  diameter,  through  the  heat 
pipe  for  an  actuation  rod.  Figure  3  illustrates  how  the  heat  pipe  is  an  integral 
part  of  the  fin. 
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Actual  Size 


Figure  2.  Heat  pipe  fin  cross-section. 


1.3  Aerodynamic  Heating  &  Cooling  Estimates 

The  phase  1  final  report,  by  Chang  (1995),  included  a  method  to  estimate  the  heat 
transfer  coefficient  for  the  airfoil  used  in  this  project.  Basically,  the  airfoil  is  treated  as  a 
flat  plate  except  for  the  leading  edge  portion.  This  front  part  of  the  fin  is  treated  as  a 
circular  cylinder  because  the  flat  plate  correlation  would  greatly  over  estimate  the 
convection  heat  transfer  coefficient,  h,  for  this  portion  of  the  fin.  These  estimated 
values  have  been  used  in  order  to  facilitate  temperature  dependent  h  estimates  for 
both  the  evaporator  and  condenser  portions  of  the  heat  pipe  for  various  operating 
temperatures. 

1.3.1  Circular  Cylinder  Heat  Transfer  Coefficient 

The  main  challenge  is  to  find  the  appropriate  diameter  of  the  circular  cylinder  used 
for  this  estimate.  Through  graphical  and  computational  analysis,  a  diameter  of  0.0127 
m  was  used  for  this  estimate.  This  circular  cylinder  correlation  is  used  for  the  first 
0.007  m  along  the  airfoil  surface  from  the  leading  edge  stagnation  point.  This 
translates  to  63°  of  the  cylinder.  The  average  h  for  the  leading  edge  is  approximated 
by  the  following  equation: 


Figure  3.  Simplified  illustration  of  heat  pipe  fin. 
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1.3.2  Flat  Plate  Heat  Transfer  Coefficient 

The  remaining  surface  of  the  air  foil,  0.007  m  to  0.116  m,  is  assumed  to  be  a  flat 


plate.  The  average  heat  transfer  coefficient,  h,  for  this  region  is  estimated  from 


hfp  =h 


0.116m 
X  0.007m 


0.0542(k) 


nO.8 

pu  ' 


1.3.3  Combined  Estimate  of  Heat  Transfer  Coefficient 
In  order  to  illustrate  how  well  the  0.0127  m  diameter  circular  cylinder  correlation 
matches  up  to  the  flat  plate  correlation  from  x=0.007  m.  Figure  4  shows  h  as  a  function 
of  surface  coordinate  distance  from  the  leading  edge  stagnation  point  for  an 
evaporator  temperature  of  900  K  . 
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Figure  4.  Heat  transfer  coefficient,  h,  profile  for  an 
evaporator  temperature  of  900  K. 


The  average  convection  heat  transfer  coefficient,  have,  is  of  course  the  weighted 

average  of  the  hcyi  and  hfp  given  above. 

i-  0.007;-  0.109  r- 

have  =  -rTTTT  ncyl + ttTTT: 


0.116 


0.116 
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2.  Heat  Pipe  Background 

2.1  Theoretical  Performance  Limitations 

There  are  many  limitations  on  the  operating  capability  of  a  heat  pipe,  but  the  major 
constraint  is  usually  the  capillary  limit.  The  entrainment  and  sonic  limits  are 
sometimes  relevant.  Chang  (1995)  gave  a  good  tutorial  on  these  limits  in  the  phase  1 
report,  and  a  detailed  description  can  be  found  in  Chi  (1976). 


2.1.1  Capillary  Limit 


-pig(Le+Lc) 


H^(Fi+Fv) 


where  rc  is  the  effective  capillary  radius,  pi  is  the  liquid  density,  g  is  the  acceleration 
due  to  gravity,  Le  and  Lc  are  the  evaporator  and  condenser  lengths.  F|  and  Fy  are 
frictional  coefficients  for  the  liquid  and  vapor  sections: 


F|  =  — ^ - 

^wPl^fg 

P  _  (fy  R6v)P'V 

2AyryPyhfg 


2.1.2  Sonic  Limit 


QS  ~  AyPyhfg 


7RTy 
2(y  +  1) 


n1/2 


2.1.3  Entrainment  Limit 


|1/2 


^ent  ~  ^v^fg 


gpv 

2>'h,s 


2.2  Steady  State  Energy  Balance 

At  steady  state,  the  energy  into  the  device  must  be  the  same  as  the  energy  taken 
out  of  the  device.  For  this  situation,  the  heat  into  the  evaporator  must  be  equal  to  the 
heat  rejected  at  the  condenser.  The  heat  into  the  evaporator,  Q,  is  given  by 

Q  =  pLehe(Tr  e  ~  Tg) 

where  p  is  the  perimeter  of  the  fin,  hg  is  the  convection  heat  transfer  coefficient  for  the 
evaporator,  Tr,e  is  the  recovery  temperature  at  the  evaporator  surface,  and  Tg  is  the 
outer  wall  temperature  for  the  evaporator.  A  similar  expression  can  be  written  for  the 
condenser  side  and  equated  to  the  evaporator  expression  yielding 

hgAg  (Tr^g  -  Tg )  =  hgAg  (Tg  -  Tr,c ) 
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After  dividing  through  by  p,  Lc  can  be  resolved  in  terms  of  Le: 


i-c  -  Lg 


"h.YT,e-Tg 


Tc-Tr,, 


oj 


The  temperature  drop  across  the  fin  wall  and  the  saturated  wick  is  calculated 
based  on  one-dimensional  conduction  of  a  steady  state  heat  load,  Q.  The  temperature 
of  the  inside  wall  of  the  heat  pipe  evaporator,  Tw,e.  is  given  by 


T 


w,e 


=  Te 


Qtp 

pLgkp 


where  tp  is  the  wail  thickness,  and  kp  is  the  thermal  conductivity  of  the  heat  pipe 
container  wall.  Similarly,  the  vapor  temperature,  Ty,  is  calculated  from 

T  —  T  Q^wick 

'  V  ~  '  w,e  „t  I, 

pLgkgff 

where  twick  is  the  wick  thickness,  and  kgff  is  the  effective  thermal  conductivity  of  the 
saturated  wick. 

Unfortunately,  the  condenser  temperature,  Tc,  is  still  an  unknown  in  these 
equations.  It  can  be  determined  by  the  following  relationship: 

^wick  , 
kgff  kp 

It  can  be  easily  seen  that  the  above  equations  are  coupled,  but  what  is  not  obvious 
is  that  many  of  the  equations  have  terms  that  are  nonlinear  functions  of  the 
temperatures  (like  the  h's  and  k's).  So  this  ends  up  to  be  a  set  of  nonlinear  coupled 
equations  that  must  be  solved  in  an  iterative  manner.  Since  it  is  desired  that  this  heat 
pipe  operates  at  as  low  a  temperature  as  possible,  solution  of  these  equations  is 
automated  using  a  Microsoft  Excel  spreadsheet  for  several  potential  working  fluids. 


2.3  Operating  Temperature  ....  Working  Fluids 

The  desired  operating  temperature  for  this  heat  pipe  is  in  the  600  -  700  K  range. 
But  from  previous  experience  this  seems  improbable  given  the  existing  working  fluids 
for  this  temperature  range.  The  following  sections  outline  the  process  of  determining 
the  best  heat  pipe  working  fluid  for  a  given  temperature  range.  This  is  done  by 
estimating  the  heat  transfer  coefficient  at  the  evaporator  and  condenser,  then 
determining  the  total  heat  input  to  the  evaporator,  and  then  estimating  the  condenser 
length.  These  are  all  functions  of  temperature. 

2.3.1  Average  Heat  Fluxes 

The  evaporator  and  condenser  heat  transfer  coefficients  have  been  approximated 
using  flat  plate  and  circular  cylinder  correlations.  The  leading  edge  of  the  air  foil,  from 
0  to  0.007  m,  is  estimated  to  be  a  circular  cylinder  of  0.0127  m  diameter;  and  the  rear 
of  the  airfoil  is  approximated  as  a  flat  plate.  The  heat  transfer  coefficient,  h,  has  a 


5 


maximum  value  at  the  leading  edge  and  steadily  decreases  with  the  chord  length  (see 
Figure  4).  The  calculated  average  heat  transfer  coefficients  are  plotted  in  Figure  5. 

Now  given  values  for  h  as  a  function  of  evaporator  and  condenser  temperatures,  Q, 
q,  and  Lc  can  be  calculated  as  functions  of  temperature.  These  quantities  are 
important  in  determining  a  feasible  working  fluid  and  operating  temperature. 


Figure  5.  Calculated  average  heat  transfer  coefficients  as  a 
function  of  temperature  for  the  evaporator  and  condenser. 

The  heat  flux  into  the  evaporator  section  of  the  heat  pipe  is  a  function  of  the  heat 
pipe  operating  temperature.  Using  the  estimates  for  h  given  above,  the  total  heat  input 
as  a  function  of  evaporator  temperature  has  been  calculated  and  is  illustrated  in 
Figure  6.  As  expected,  the  heat  input  to  the  evaporator  decreases  as  the  evaporator 
temperature  increases. 


Figure  6.  Thermal  input  power,  Q,  as  a  function  of 
evaporator  temperature. 
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2.3.2  Evaporator  Average  and  Maximum  Heat  Fluxes 

The  average  and  maximum  (leading  edge)  heat  fluxes  have  been  calculated  and 
plotted  as  a  function  of  evaporator  temperature  in  Figure  7.  As  can  be  seen  from 
Figure  7,  the  average  heat  fluxes  are  veiy  high  below  800  K  and  are  likely  to  be 
beyond  the  capability  of  most  known  working  fluids  for  that  temperature  range.  The 
leading  edge  heat  fluxes  are  typically  about  twice  the  average  fluxes,  making  the 
situation  even  worse. 


Figure  7.  Leading  edge  and  average  heat  fluxes  into  the 
evaporator  as  a  function  of  evaporator  temperature. 


2.3.3  Condenser  Length 

Even  if  a  working  fluid  can  be  found  to  operate  in  the  desired  temperature  range 
and  handle  the  heat  fluxes,  there  may  still  be  a  practical  limit  on  the  length  of  the 
condenser  given  that  this  is  a  short  fin.  Condenser  lengths  for  various  evaporator 
temperatures  have  been  calculated  and  plotted  in  Figure  8.  It  can  be  seen  from  this 
figure  that  at  600  K,  the  condenser  length  would  be  about  0.2  m.  Given  that  the 
evaporator  length  is  fixed  at  about  0.5  m,  it  may  be  feasible  to  operate  at  600K  if  a 
suitable  working  fluid  can  be  found. 


Figure  8.  Calculated  condenser  lengths  as  a 
function  of  evaporator  temperature. 
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2.3.4  Working  Fluids 

Now  that  it  has  been  shown  from  an  energy  balance  that  the  fin  is  long  enough  to 
accommodate  a  600  -  700  K  heat  pipe,  it  must  be  determined  if  there  is  a  working  fluid 
that  can  do  the  job.  Figure  9  shows  most  of  the  useful  heat  pipe  working  fluids  and 
their  useful  operating  temperature  ranges.  It  can  be  seen  from  this  figure  that  Hg,  Cs, 
and  H2O  are  the  only  possibilities  at  600  -  700  K. 


I - 1 _ I _ I _ LJ _ ■  I  "I  I  I  I  I  I  I  I  I  I 

-440  -400  -300  -200  -100  0  100  500  1000  10000 

TEMPERATURE,  “F 

Figure  9.  Heat  pipe  working  fluids  and  operating  temperatures  (Chi,  1976). 

Now,  given  that  these  three  working  fluids  can  operate  in  this  temperature  range,  it 
must  be  determined  if  any  of  them  can  handle  the  required  heat  loads.  One  indicator 
of  the  heat  transport  capability  of  a  fluid  is  the  liquid  transport  factor,  N|,  defined  as  the 
surface  tension  times  the  latent  heat  divided  by  kinematic  viscosity: 


Figure  10  shows  N|  as  a  function  of  temperature  for  water,  Hg,  and  Cs.  Water  is 
clearly  an  unacceptable  working  fluid  near  or  above  600  K.  This  fact  is  not  surprising 
because  water  has  a  critical  temperature  of  about  643  K.  It  can  also  be  seen  that  Cs 
has  a  very  low  N|  in  general,  making  it  an  unacceptable  working  fluid  for  this  high  heat 
flux  application.  It  does  appear  that  Hg  has  a  N|  that  could  be  acceptable. 
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The  next  indicator  to  look  at  is  the  wicking  height  factor,  H.  This  factor  is  defined  as 
the  ratio  between  the  surface  tension  and  the  density  times  gravitation  acceleration. 
This  term  gives  an  indication  of  how  great  a  pressure  head  can  be  supported  in  the 
wick  of  the  heat  pipe.  This  term  is  important  for  this  application  given  that  the  working 
fluid  must  wick  up  against  gravity  in  the  bottom  fin.  This  term  is  given  by 


Pig 


Figure  10.  Liquid  transport  factors  as  a  function  of  fiuid 
temperature  (Data  from  B  &  K  Volume  il). 


Figure  1 1  plots  H  as  a  function  of  fluid  temperature  for  water,  Hg,  and  Cs.  Again,  it 
is  clear  from  this  figure  that  water  is  unacceptable  anywhere  near  600  K. 
Unfortunately,  Cs  and  Hg  also  have  poor  values  of  H  for  the  required  wicking  distance 
of  this  application. 


Figure  11.  Wicking  height  factors  as  a  function  of  fluid 
temperature  (Data  from  B  &  K  Volume  II). 


2.3.5  Selection  of  Working  Fluid 

From  the  above  discussion,  it  may  appear  that  Cs  and  Hg  are  not  acceptable 
however,  calculations  must  be  done  to  be  certain.  The  calculations  were  run  for  both 
working  fluids,  and  they  all  are  capillary  limited  at  all  temperatures  calculated. 
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After  these  calculations  were  accomplished,  there  was  only  the  remote  possibility 
that  one  of  Mainstream's  new  working  fluids  could  do  the  job.  Mainstream,  a  Florida 
based  thermosciences  company,  was  given  the  heat  flux  requirements  and  other 
pertinent  physical  requirements  to  evaluate  if  one  of  these  new  fluids  would  be 
suitable  for  this  application.  Mainstream  estimates  that  none  of  these  new  fluids  is 
close  to  having  the  required  transport  and  wicking  ability  for  this  application. 

This  left  no  choice  but  to  use  either  Na  or  K  as  the  working  fluid;  and  this  fact 
imposes  that  the  operating  temperature  must  be  above  800  K.  At  900  K,  Na  has  a  N|  of 

2.2x1 0'*2  W/m2  and  an  H  of  1.76x10’^  m^,  which  should  be  adequate  for  this 
application. 
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3.  Could  Insulation  Help? 

Don  Thomas  of  the  US  Army  Missile  Command  raised  an  interesting  question:  If 
enough  insulation  were  placed  between  the  heat  source  and  the  heat  pipe  working 
fluid,  could  the  heat  pipe  operate  at  a  temperature  suitable  to  use  water  as  the  working 
fluid?  This  is  certainly  an  interesting  question  and  worth  exploring  further. 

To  facilitate  answering  this  question,  consider  the  simplified  rnodel  of  this 
hypothetical  conduction  path  illustrated  in  Figure  12.  Accounted  for  here  are  the 
conduction  paths  for  the  stainless  steel  fin  wall  (kfjn,  Xfjn),  insulator  material  (kins,  Xins), 
heat  pipe  wall  (khp,  Xhp),  and  heat  pipe  wick  (kwick,  Xwick)-  Neglected  are  all  of  the 
contact  resistances  between  the  material  layers.  For  a  steady  state  situation,  the  heat 
flow  (q)  through  each  layer  in  Figure  12  will  be  the  same: 

^  =  ^(Ts-T,)  =  |^(T,-T2)  =  ^(T2-T3)  =  ^(T3-T„) 

'  ^  Xfjn  X|ns  Xhp  x^ick 

where  Ts  is  the  outer  surface  temperature  of  the  fin,  and  Ty  is  the  vapor  temperature 
of  the  heat  pipe  as  before.  The  intermediate  temperatures  can  easily  be  eliminated: 

n/  (Ts-Ty)  (Ts-T.) 

Xfin  I  Xjns  I  I  Xwick  1^ 
l^fin  *^ins  l^hp  kv(,|ck 

The  heat  flux  input  is  determined  by  the  free  stream  temperature  (Too)  and  the  free 
stream  convection  coefficient  (hoo): 

This  equation  can  now  be  solved  for  Ts' 


T  _i_ 

r  1  ^ 

1  +  heoD  j 

1  oo  + 

j^1*fheoD  j 

Now  the  fin  surface  temperature  Ts  can  be  easily  solved  as  a  function  of  insulator 
thickness,  Xjns,  if  all  other  parameters  are  held  constant.  For  this  calculation,  the  vapor 


Xfin 

Xjns 


^hp 

Xwicl4 


fin  skin,  kfjp 
insulator,  kjng 
heat  pipe  wall,  khp 
wick,  k^fck 


Figure  12.  Simplified  conduction  path 
including  an  insulation  layer. 
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temperature  (Tv)  will  be  fixed  at  473  K,  close  to  the  upper  practical  operating  limit  for 
water.  Values  of  1200  W/m^-K  for  hoo  and  1213  K  for  To©  are  used  as  in  the  previous 
heat  pipe  calculations.  The  other  parameters  used  for  this  calculation  are  shown  in 
Table  1. 


Table  1.  Parameters  Used  for  Calculation. 


X  (m) 

k  (W/m-K) 

material 

fin  skin 

0.0023 

22 

stainless  steel 

insulator 

variable 

0.058 

Fiberfrax  Ceramic  Fiber 
Paper  (HSA) 

heat  pipe  wall 

0.0023 

375 

copper 

wick 

0.00023 

1.22 

copper/water 

Figure  13  shows  both  the  fin  skin  temperature  and  the  input  heat  flux  to  the 
evaporator  as  a  function  of  the  insulator  thickness  for  a  heat  pipe  vapor  temperature 
fixed  at  473  k.  In  this  simplified  calculation,  the  fin  surface  temperature  increases 
rapidly  as  with  increases  in  insulator  thickness.  As  the  fin  skin  temperature  increases, 
the  heat  flux  into  fin  rapidly  decreases  because  this  is  driven  by  the  temperature 
difference  between  the  fin  and  the  free  stream  temperature.  There  would  of  course  be 
a  practical  upper  limit  to  the  fin  skin  temperature;  and  for  stainless  steel  this  would  be 
around  1000  K.  It  can  be  seen  in  Figure  13  that  the  insulator  thickness  corresponding 
to  a  1000  K  fin  temperature  is  just  slightly  above  0.0001  m,  predicting  a  heat  flux  of 


Figure  13.  Fin  skin  temperature  and  thermal  input  heat 
flux  as  a  function  of  insulator  thickness. 


12 


about  25  W/cm2. 

These  calculations  are  clearly  based  on  assumptions  that  overly  simplify  the  actual 
physical  situation,  but  the  results  suggest  that  the  use  of  a  water  heat  pipe  in  this  type 
of  situation  Is  feasible.  A  copper/water  heat  pipe  is  certainly  feasible  near  25  W/cm^  or 
below,  and  the  insulator  thickness  required  to  get  this  heat  flux  is  only  minimal.  This 
idea  certainly  warrants  further  more  thorough  exploration. 
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4.  Heat  Pipe  Design  Specifications 

Thermacore  Inc.,  a  small  business  in  Lancaster  PA  specializing  in  heat  pipe 
technology,  was  given  the  following  specifications  to  fabricate  a  stainless  steel  -  Na 
heat  pipe  to  operate  above  800  K  with  the  thermal  input  and  heat  fluxes  indicated  in 
Figure  14.  Table  2  contains  relevant  heat  pipe  parameters. 


Figure  14.  Heat  flux  and  thermal  input  information. 


Table  2.  Heat  Pipe  -  Fin  Parameters 


Item 

Specification 

Operatinq  Temp 

800  -  900  K 

Wall  Material 

304  Stainless  Steel  (SS) 

Working  Fluid 

Sodium  (Na) 

Wick 

200  mesh  304  SS  screen 

Configuration  Shape 

NACA0015  (c=0.1143  m) 

Wall  Thickness 

0.0023  m 

Evaporator  Length 

0.050  m 

Condenser  Length 

0.055  m 

Input  Power 

3,186  W  @Te  =  850K 

Evaporator  Heater 

Radiant  removable 

Condenser  Cooler 

heat  choked  water 

Total  Fin  Length 

0.254  m 

Acceleration  Field 

Short  duration  3g 

Thermacore  built  and  delivered  two  proof-of-concept  heat  pipe  fins.  Appendix  A 
contains  a  fabrication  summery  written  by  Thermacore,  and  Appendix  B  has  the 
fabrication  and  assembly  drawing  that  Thermacore  used.  Both  heat  pipes  failed 
shortly  after  steady-state  testing  due  to  a  faulty  weld. 


14 


5.  Transient  Analysis 

An  analytical  model  similar  to  those  developed  by  Silverstein  (1971),  Camarda 
(1977),  and  Chang  (1996)  was  used  to  predict  the  transient  operation  from  frozen 
start-up  of  a  liquid  metal  heat  pipe.  The  analytical  model  consists  of  lumped 
capacitance  models  representing  the  five  possible  physical  stages  encountered 
during  the  start-up  of  a  liquid  metal  heat  pipe  from  a  frozen  state.  Figure  15  depicts  the 
various  stages  of  start-up  described  in  detail  below. 


Qe 


:e - ►! 

-Zn - ► 

- Zn+1 - ^ 


Figure  15.  Control  volumes  used  for  transient  analysis. 


5.1  Stage  1:  Frozen  Working  Substance 

This  stage  predicts  the  transient  temperature  Increase  of  the  evaporator  from  heat 
input  to  the  evaporator.  The  assumptions  for  this  stage  are:  a)  the  evaporator  can  be 
treated  as  one  lumped  capacitance  body;  b)  only  the  evaporator  changes  in 
temperature  with  heat  input  to  the  evaporator;  c)  the  working  substance  remains  in  a 
frozen  state  during  this  stage.  Therefore,  for  T^p  <  T^eit.  the  new  temperature  of  the  heat 
pipe  after  a  time  step.  At,  is  predicted  by 


*T»n+l _ rpH 

^hp  ~  ^hp 


QeAt 

CjZe 


(  1  ) 


where  Ci  =  (pCp)^  Ap  -h  (1  -  e)(pCp)^  Av^  +  £(pCp)^  A^ . 


5.2  Stage  2:  Evaporator  Thaw 

This  stage  predicts  the  time  required  to  thaw,  the  working  substance  in  the 
evaporator.  The  assumptions  for  this  stage  are:  a)  the  evaporator  can  be  treated  as 
one  lumped  capacitance  body;  b)  only  the  evaporator  is  heated;  c)  there  is  no  change 
in  temperature.  Therefore,  for  T^p  =  T,„en,  the  time  to  melt  the  evaporator  section  is 
given  by 


Atmelt  ”■  ^Pl Aw^e  p. 

Ve 


(2) 


15 


5.3  Stage  3:  Evaporator  Heating  to  Transition  Temperature 

This  stage  predicts  the  transient  temperature  increase  of  the  thawed  evaporator  up 
to  the  transition  temperature,  T„  of  the  working  fluid,  predicted  by 


=  4990 


7C 

R 


M'v 

vPvly 


(  3  ) 


Due  to  the  temperature-dependent  properties  of  the  parameters  used  in  equation 
(3),  the  transition  temperature  must  be  calculated  using  an  iterative  process.  The 
assumptions  for  this  stage  are:  a)  the  evaporator  can  be  treated  as  one  lumped 
capacitance  body;  b)  only  the  evaporator  changes  in  temperature  with  heat  input  to 
the  evaporator;  c)  the  vapor  space  is  a  free  molecular  region,  thus  heat  transport  out 
of  control  volume  a-b  in  Figure  15  is  neglected.  Therefore,  for  T^^g,,  <  T(,p  <  T„  the  new 
temperature  of  the  heat  pipe  after  a  time  step,  At,  is  predicted  by 

tH+I  _  I  /  >1  \ 

where  Cj  =  (pCp)  Ap  +(1  -e)(pCp)^A„ +e(pCp)|A„. 


5.4  Stage  4:  Meit  Front  Propagation 

In  this  stage,  heat  pipe  operation  begins  as  the  continuum  front  moves  down  the 
length  of  the  heat  pipe.  The  assumptions  for  this  stage  are:  a)  the  continuum  region, 
length  z",  can  be  treated  as  one  lumped  capacitance  body;  b)  the  continuum  region  is 
established  initially  over  the  entire  evaporator  region,  Zg,  as  soon  as  the  transition 
temperature  is  reached;  c)  T,  <  T^p  and  Zg  <  z"  <  Zg+z^.  An  energy  balance  of  control 
volume  b-c  in  Fig.  15  gives 


Qf  =  Qs  -  •>cp(z"  -  z.)(tJp  -  T..)  -  ■— - 
Similarly,  an  energy  balance  of  control  volume  a-c  gives 

Qe  =  Qf  -h  hcp(z"  -  Ze)(Tfp  -  T„)  +  -  T^,) 


(  5  ) 

(6) 


Finally,  an  energy  balance  on  control  volume  c-d  in  gives  the  next  continuum  front 
position  after  time  step.  At,  given  by 

- r -  (7) 

”  Xjo)  +  C2|T  —  j  *+•  psA^hg] 

Equations  (5)  and  (6)  are  coupled  and  nonlinear  due  to  temperature-dependent 
properties,  so  an  iterative  solution  process  must  be  used  to  solve  for  Q,  and  using 
the  current  continuum  front  position,  z".  Equation  (7)  is  then  solved  once  Q,  and 
have  been  determined.  This  stage  is  continued  until  the  continuum  front  reaches  a 
steady-state  position,  or  the  continuum  front  reaches  the  end  of  the  condenser  (stage 
5). 
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5.  Stage  5:  Condenser-Limited  Heat  Pipe  Operation 

In  this  stage,  the  entire  heat  pipe  is  active  and  the  transient  temperature  increase  is 
predicted  until  steady-state  temperature  is  reached.  The  assumptions  for  this  stage 
are:  a)  the  entire  heat  pipe  can  be  treated  as  one  lumped  capacitance  body;  b)  the 
continuum  region  is  established  over  the  entire  heat  pipe  length,  z"  =  Zg+z^.  The  new 
temperature  of  the  heat  pipe  after  a  time  step.  At,  is  predicted  by 

rpn+l  -T-n  .  (Qe~Qc)^*-  f  8  ) 


5.6  startup  Predictions 

The  procedures  and  stages  for  predicting  the  transient  startup  of  a  liquid  metal  heat 
pipe  outlined  in  equations  (1)  -  (8)  were  implemented  using  a  Mathematica  2.2 
program.  A  program  listing  and  sample  calculation  is  provided  in  Appendix  C. 
Mathematica  was  chosen  to  implement  this  model  due  to  its  built-in  iteration 
capabilities  and  solvers  for  systems  of  coupled  (even  non-linear)  equations. 

Figure  16  is  the  transient  prediction  for  the  front  heat  pipe.  The  top  portion  of  the 
figure  is  a  prediction  of  the  evaporator  temperature  until  the  transition  temperature 


Figure  16.  Transient  prediction  for  front  pipe. 
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(832  K)  is  reached,  then  it  is  the  temperature  of  the  continuum  region.  As  shown  here, 
the  model  has  assumed  that  the  heat  pipe  will  operate  at  the  transition  temperature 
unless  the  condenser  length  is  too  short.  The  bottom  portion  of  Figure  16  is  the 
predicted  continuum  front  propagation  after  start  up.  It  can  be  seen  here  that  there 
should  be  no  continuum  front  until  the  transition  temperature  is  reached,  then  the  front 
begins  to  propagate  past  the  evaporator  section.  The  steady-state  time  predicted  to  be 
205  seconds  by  determining  where  this  curve  reaches  its  asymptotic  maximum.  Note 
that  since  this  curve  does  appear  horizontal  at  steady-state,  the  condenser  length  is 
long  enough.  If  the  condenser  length  is  too  short,  the  continuum  front  verses  time 
curve  would  stop  short  of  horizontal  as  the  end  of  the  condenser  is  reached,  and  the 
temperature  would  increase  above  the  transition  temperature. 

A  similar  prediction  can  be  seen  in  Figure  17  for  the  rear  heat  pipe.  The  rear  heat 
pipe  is  predicted  to  reach  a  steady  state  operating  temperature  of  822  K  after  212 
seconds. 


Figure  17.  Transient  prediction  for  rear  pipe. 
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6.  Results  and  Conclusions 

Two  proof-of-concept  heat  pipe  fins,  one  shown  In  Figure  18,  were  fabricated  by 
Themiacore  and  verified  for  steady-state  operation.  Appendix  A  contains  a  fabrication 
summery  written  by  Thermacore,  and  Appendix  B  has  the  fabrication  and  assembly 
drawings  that  Themiacore  used.  Both  heat  pipes  failed  shortly  after  steady-state 
testing  due  to  a  faulty  weld.  As  a  result,  tiie  transient  predictions  cannot  be  validated 
and  the  use  of  a  liquid  metal  heat  pipe  for  this  translent-in-nature  application  cannot 
be  fully  proven.  It  should  be  noted  that  this  failure  was  due  to  a  fabrication  problem 
and  not  a  design  problem. 


Figure  18.  Proof-of-concept  heat  pipe  fin. 


6.1  Steady-State  Results 

The  experimental  results  of  the  steady  state  testing  are  given  in  detail  in  Appendix 
A,  but  some  of  the  more  pertinent  observations  follow: 


Case 

1 

2 

3 

4 

5 

6 

average 


Table  3.  Experimental  steady-state  results. 


Predicted  SS 
Temp  (K) 
820 
820 
820 
820 
820 
820 
820 


Actual  Front 
HP  Temp  (K) 
839 
852 
820 
819 
832 
793 

825.83 


Actual  Rear 
HP  Temp  (K) 
854 
881 
831 
829 
829 
825 
841.5 
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It  can  readily  be  seen  from  Table  3  that  the  predicted  steady-state  temperature  is 
within  1  %  of  the  average  experimental  value  for  the  vapor  space  thermocouple  on  the 
front  heat  pipe,  and  within  3%  of  the  average  vapor  space  value  for  the  rear  heat  pipe. 
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7.  Cost  Summary 

The  following  is  a  breakdown  of  how  the  funds  received  by  MIPR  from  the  Army 
Missile  Command  were  used: 


Labor  hours 

$43K 

Heat  Pipe  Fabrication 

$39K 

Test  equipment 

$  8K 

$90K 

\ 
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3.0  TECHNICAL  APPROACH 


This  section  of  the  report  describes  the  technical  approach  used  to  meet  the  heat  pipe  fm 
specifications  listed  in  the  solicitation,  as  summarized  in  Appendix  A.  Each  specification  is  addressed 

in  the  following  discussion. 

1.  Operating  temperature  between  800K  and  900K. 

Thermacore  has  completed  a  preliminaiy  analysis  of  the  missile  fin  heat  pipe,  based  on  the 
information  given  in  Appendix  A.  Due  to  the  size  of  the  control  hole,  the  fin  was  divided  into  two 
separately  operating  heat  pipes:  one  in  firont  of  the  control  hole  and  one  behind  the  control  hole.  A 
calculation  was  made  to  estimate  the  sonic  limit  for  the  fi-ont  heat  pipe,  which  will  be  subjected  to  the 
larger  axial  heat  flux.  It  was  found  that  the  front  heat  pipe  must  rise  to  a  temperature  of  about  900K 
for  conservative  operation  of  the  front  heat  pipe  below  the  sonic  limit  to  occur.  Thus,  operation  of  this 
device  under  simulated  operating  heating  conditions  probably  will  not  allow  normal  operation  of  the 
front  heat  pipe  at  800K,  although  the  required  test  conditions  will  be  at  lower  leading  edge  heat  fluxes 
than  actual  operating  heat  fluxes,  so  this  limit  may  not  be  observed  in  testing. 

2.  Wall  material. 

Thermacore  assembled  the  missile  fin  proof-of-concept  device  from  0.090"  thick  304L 
stainless  steel  sheet,  joined  primarily  by  welding.  Sodium  heat  pipes  using  welded  formed  sheet 
construction  were  designed  and  built  by  Thermacore  under  Air  Force  Contract  F33615-86-C-3227,  as 
described  below  under  References  and  Proof  of  Performance. 


3.  Working  fluid. 

Sodium  was  the  selected  working  fluid.  It  has  been  successfully  used  in  several  leading  edge 
heat  pipes  previously  built  by  Thermacore.  In  addition  to  the  work  described  in  Item  2  above, 
Thermacore's  large  solar  receiver  heat  pipe  program  is  based  on  sodium  heat  pipe  technology. 
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4.  Wick  Material. 

A  porous  metal  nickel  (-50  +  100)  wick  structure  rather  than  a  stainless  steel  screen  wick  was 
selected  for  this  application.  Theimacore  has  had  good  success  with  this  type  of  wick  structure  for  our 
Solar  Applications  products,  i.e.  sodium/superalloy  heat  pipe  solar  receivers.  Contracts  listed  in  the 
reference  list  provide  several  examples.  Absorbed  heat  fluxes  of  up  to  lOOW/cm^  are  regularly 
achieved  with  sodium  working  fluid  using  such  wicks.  In  addition,  sintered  porous  wicks  with  a 
distribution  of  pores  have  been  analytically  and  experimentally  shown  to  be  boiling  tolerant,  in  contrast 
to  single  pore  size  wicks  such  as  layered  screens.  This  approach  will  thus  provide  a  favorable 
performance  margin  for  the  missile  fin  application.  Preliminary  analysis  has  shown  that  a  twofold 
safety  margin  can  be  realized  using  this  design.  Wick  properties  for  the  recommended  wick  structure 
are  summarized  below. 

permeability  =  3.5E-llm^ 

pore  radius  =  4.0  E-5  m 

thermal  conductivity  =  60  W/m/K 

wick  thickness  =  1.25  mm 

5.  Configuration  shape  NACA  0015  (c=0.1 143  m). 

Theimacore  formed  the  configuration  shape  using  techniques  previously  used  successfully  for 
related  leading  edge,  turboshaft  engine,  and  solar  receiver  geometries.  The  leading  edge  radius  was 
shaped  around  a  forming  mandrel.  The  remainder  of  the  fin  shape  was  set  by  joining  the  cover  sheet 
to  machined  end  plates  and  internal  stiffeners  which  were  machined  to  the  specifications  listed  in 
Appendix  A  of  the  solicitation.  Theimacore  also  has  an  active  program  that  uses  the  same  approach 
for  forming  turboshaft  engine  vane  heat  pipes  which  have  a  similar  shape  to  those  required  for  the 
missile  fin.  This  work  is  being  done  on  US  Army  Contract  DAAJ02-94-C0023,  and  also  uses  the  same 
forming  technologies  proposed  for  use  on  the  subject  proposed  work. 

Instrumentation  for  testing  the  heat  pipe  included  thermocouples  to  measure  the  vapor 
temperature  of  each  heat  pipe  and  to  perform  calorimetry  on  the  test  article.  An  instrumentation  and 
assembly  drawing  was  prepared  and  submitted  to  Wright  Laboratory  for  review  and  approval  prior  to 
fabrication  of  the  hardware. 
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6.  Wall  thickness  0.0023  m. 

The  device  was  fonned  from  sheet  0.0023  m  (0.090")in  thickness,  as  described  in  Item  2  above. 

7.  Evaporator  length  0.050  m. 

Sizing  calculations  were  performed  using  this  evaporator  length,  which  showed  that  the 
recommended  approach  will  successfully  handle  the  heat  flux,  power  levels,  and  operating 
temperatures  specified  herein\ 

8.  Condenser  length  0.072  m. 

The  same  analysis  described  in  Item  7  above  showed  that  the  front  heat  pipe  and  possibly  the 
rear  heat  pipe  will  require  a  condenser  length  longer  than  5.5  cm  to  function  as  a  normal  heat  pipe  at 
prescribed  conditions  and  convect  heat  from  the  condenser  surfaces.  Lengths  of  the  heat  pipe 
condensers  were  specified  from  detailed  analyses  performed  during  the  final  design,  as  shown  in 
Appendix  C. 

9.  Input  power  3 . 1 86  kW  at  850K. 

Pieliminaiy  analysis  showed  that  the  selected  design  can  handle  the  required  input  power.  A 
radiant  heat  source  was  used,  which  is  capable  of  supplying  the  required  input  power,  as  described 
under  Item  10. 

1 0.  Evaporator  heat  radiant  removable  type. 

Thermacore  used  quartz  lamp  heaters  for  testing  the  heat  pipes.  Thermacore  has  developed 
considerable  experience  using  this  approach  for  the  solar  receiver  heat  pipe  programs  described  below 
under  related  work.  A  procedure  for  operating  and  testing  the  heat  pipe  was  supplied  with  the 
hardware  at  the  time  of  delivery. 

11.  Condenser  cooler  affixed  heat  choked  water-cooled  type. 

Thermacore  built  and  delivered  a  water-cooled  clamp-on  type  condenser/calorimeter  which  was 
supplied  with  the  heat  pipe  at  time  of  delivery.  Standoffs  were  used  on  the  condenser  so  that  the  water 
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in  the  condenser  does  not  boil  and  so  that  the  heat  pipe  is  able  to  operate  normally  without  freezing 
the  woridng  fluid  or  "shocking"  die  heat  pipe.  This  latter  non-operational  condition  occurs  when  the 
condenser  design  is  too  efficient,  and  prevents  the  heat  pipe  condenser  from  operating  at  subsonic 
conditions  during  startup.  Previous  experience  in  design,  fabrication,  and  testing  of  similar  devices 
was  used  to  produce  a  condenser  capable  of  allowing  the  heat  pipe  to  operate  in  conditions  simulating 
normal  operation. 

12.  Total  fin  length  0.254  m. 

The  total  fin  length  and  geometry  was  built  based  on  the  specifications  listed  in  Appendix  A, 
as  attached.  Sheet  metal  extensions  were  added  to  the  heat  pipes  to  form  the  full  length  of  the  fin. 

1 3 .  Acceleration  field  3g  for  short  duration. 

A  preliminary  analysis  was  performed  which  showed  that,  during  the  20  second  acceleration 
at  2g,  the  temperature  rise  of  the  heat  pipe  will  not  be  enough  to  raise  the  heat  pipe  to  operating 
temperature.  Thus,  acceleration  effects  on  the  heat  pipe  should  be  minimal  except  during  the  final 
seconds  of  flight.  The  latter  condition  also  occurs  over  such  a  short  time  duration  that,  based  on 
preliminary  calculations,  the  heat  pipe  will  not  be  adversely  affected.  These  analyses  are  provided  in 
Appendix  C. 

14.  Evaporator  location  along  portion  of  fin  adjacent  to  missile  body. 

The  recommended  heat  pipe  design  is  based  on  this  specification. 

4.0  DESIGN  AND  FABRICATION  OF  THE  MISSILE  FIN  HEAT  PIPES 

Thetmacore  completed  a  series  of  performance  calculations  to  verify  the  specifications  supplied 
by  the  Air  Force  and  to  size  the  heat  pipe  wick  structure.  Table  1  summarizes  the  as-built  heat  pipe 
specifications;  the  supporting  calculations  used  to  develop  this  design  are  given  in  Appendix  C. 
Figure  4  shows  the  assembly  drawing  developed  from  the  selected  design.  The  design  consists  of  two 
separate  heat  pipe  vapor  chambers.  One  heat  pipe  is  located  in  front  of  the  control  hole;  the  second 
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TABLE  1.  -  Missile  Fin  Heat  Pipe  Design  Specifications 


Item 

Specification 

Operating  Temperature 

800  -  900  K 

Wall  material 

304  Stainless  Steel  (SS) 

Working  Fluid 

Sodium  (Na) 

Wick 

-50  +100  Sintered  powder  nickel 

Configuration  Shape 

NACA  0015  (c=0.1 143  m) 

Wall  Thickness 

0.0023  m 

Evaporator  Length 

0.050  m 

Condenser  Length 

0.072  m 

Input  Power 

3,186  W@Te  =  850  K 

Evaporator  Heater 

Radiant  removable 

Condenser  Cooler 

Heat  choked  water 

Total  Fin  Length 

0.254  m 

Acceleration  Field 

Short  duration  3g 

Evaporator  Location 

Innermost  edge  of  fin 

heat  pipe  is  located  behind  the  control  hole.  Each  heat  pipe  is  lined  with  a  0.040"  (1mm)  thick  layer 
of  -50  +100  sintered  nickel  powder  that  serves  as  a  heat  pipe  wick  structure.  The  evaporator  section 
of  each  heat  pipe  is  located  adjacent  to  the  missile  fin  body;  the  heat  absorbed  by  the  heat  pipe  from 
the  exhaust  plume  is  spread  along  the  heat  pipe  and  rejected  to  air  by  convection  and  radiation.  A 
welded  extension  is  attached  to  the  far  end  of  the  fin  'wfrich  is  not  a  heat  pipe,  but  is  simply  part  of  the 
control  fin.  The  fill  tube  ends  of  the  heat  pipes  protrude  into  the  hollow  inside  of  the  extension. 

A  fabrication  procedure  was  developed  to  assemble  the  heat  pipe  missile  fin.  Figure  5  lists  the 
main  assembly  steps  used  in  the  procedure.  Each  heat  pipe  was  assembled  with  several  stainless  steel 
tubes  inserted  into  their  vapor  spaces  to  assist  in  testing.  Each  heat  pipe  in  each  of  two  fin  test  articles 
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Figure  4.  Air  Foil  H^at  Pipe  -  Assembly 


was  processed  with  sodium  working  fluid,  and  operated  at  temperatures  of  up  to  800°C  to  vent  any 
non-condensible  gas  that  may  be  present  from  processing.  All  heat  pipes  operated  successfully  during 
these  preliminary  tests.  Figure  6  shows  two  views  of  the  assembled  heat  pipes,  and  illustrates  the  tack- 
welded  screen  pieces  used  to  attach  some  of  the  thermocouples  to  the  outer  surfaces  of  the  heat  pipes. 
The  remaining  thermocouples  were  attached  by  sliding  them  into  the  tube  inserts  installed  into  the  end 
caps  of  the  heat  pipes. 


5.0  HEAT  PIPE  TESTING 

Each  test  article  was  operated  at  temperatures  up  to  600 "C  using  a  pair  of  radiant  quartz  lamps 
for  a  heat  source  and  a  pair  of  clamped  water  calorimeters  for  a  heat  sink.  Two  views  of  the  test  setup 
are  shown  in  Figure  7.  Figure  8  shows  data  collected  for  each  heat  pipe  during  the  first  operating 
thermal  cycle  from  cold  startup,  including  a  diagram  showing  the  thermocouple  locations. 

During  the  first  thermal  cycle  after  the  hot  processing  step,  each  heat  pipe  developed  a  similar 
failure  in  the  form  of  a  wall  leak  in  the  evaporator  endcap  next  to  the  control  hole.  The  failure 
manifested  itself  as  a  non-condensible  gas  leak  into  the  front  heat  pipe,  \^hich  revealed  itself  as  a  cold 
end  in  the  condenser.  The  leak  point  was  determined  after  cooldown  to  room  temperature  by  visual 
inspection  of  the  test  articles;  an  alkali  metal  oxide  spot  formed  on  the  surface  of  the  fin  at  the  leak 
locations.  The  first  heat  pipe  was  operated  at  power  levels  up  to  1025  watts;  the  second  heat  pipe  was 
operated  at  power  levels  of  up  to  297  watts. 

The  reason  for  the  failure  is  believed  to  be  caused  by  a  combination  of  improper  wall  material, 
design  specifications,  and  fabrication  procedures.  All  previous  thin-walled  sodium  heat  pipe  wing 
leading  edges  that  were  successfully  developed  by  Thermacore  used  superalloys  rather  than  304 
stainless  steel  for  the  envelope  material,  and  all  hardware  used  a  minimum  wall  thickness  of  at  least 
0.035".  This  hardware  includes  the  Boeing  wing  leading  edges  and  the  test  articles  built  for  Wright 
Laboratory  on  the  cited  SBIR  program.  Stainless  steel  is  less  favorable  for  applications  where  vacuum 
leaktight  thin  walls  are  required,  since  it  is  more  susceptible  to  crystal  grain  growth  and  grain  cracking 
than  are  superalloys.  Furthermore,  some  of  the  procedures  used  to  assemble  the  heat  pipes  could  be 
modified 
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FABRICATION  PROCEDURE  FOR  MISSILE  FIN 


1.  CLEAN  ALL  PARTS  IN  ACETONE. 

2.  INSERT  CENTERLESS  GROUND  TUBE  IT  SHELL  AND  TACK  WELD  END  PLATES. 

3.  BRAZE  TUBE  WITH  NICROBRAZE  30  AT  1 150  oC  FOR  5  MINUTES 

4.  STOP-OFF  WELD  JOINTS. 

5.  INSPECT  BRAZE  JOINT  REBRAZE  IF  NECESSARY. 

6.  WELD  THERMOCOUPLE  WELLS  AND  LEAK  CHECK. 

7.  WELD  THERMOCOUPLE  WELLS  TO  END  PLATES. 

8.  FABRICATE  SOCK  MANDRELS  WITH  1 00  MESH  304  SS  SREEEN. 

9.  TACK  WELD  EVAPORATOR  END  PLATE  TO  AIR  FOIL. 

10.  INSERT  SOCK  MANDRELS. 

11.  FILL  WITH  -50  +100  NICKEL  POWDER  AND  VIBRATE  POWDER. 

12.  DRY  HYDROGEN  SINTER  AT  1 1 00  oC  FOR  2  HOURS. 

13.  ELECTRON  BEAM  WELD  END  PLATES  AND  TAIL  SECTION  AND  LEAK  CHECK. 

14.  BRAZE  CALORIMETER  WITH  NICROBRAZE  LM. 

15.  VACUUM  FIRE  MISSILE  FIN  FOR  1  HOUR  AT  900  oC. 

16.  CHARGE  HEAT  PIPES  WITH  SODIUM. 

17.  PROCESS  HEAT  PIPES  AT  SAME  TIME. 

18.  SEND  HEAT  PIPE  TO  ELECTRON  BEAM  WELDER  TO  HAVE  NON-HEAT  PIPE  FOIL 
WELDED  TO  HEAT  PIPE. 

19.  SETUP  LAMPS  AND  CALORIMETER  FOR  THERMAL  TEST. 

20.  ATTACH  THERMOCOUPLES  TO  THE  HEAT  PIPE. 

21.  TEST  HEAT  PIPE  ON  LAMPS  WITH  GAS-GAP  CALORIMETER. 

1 .  Turn  on  water  to  lamps  ( 2  GPM ) 

2.  T um  on  air  to  lamps  ( 1 7  CFM ) 

3.  T um  on  water  to  calorimeter  ( 250  cc’s ) 

4.  Turn  on  controller  to  1 .1  on  10  turn  pot,  let  lamps  start,  let  heat  pipe  warm  up.  Let  stabilize 
and  turn  power  supply  up  in  0.1  increments, 

5.  Let  heat  pipe  temperature  stabilize  before  turning  up  power  supply. 

1.  TAKE  PHOTOGRAPHS. 


Figure  5.  Fabrication  Procedure  for  Missile  Fin 
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Figure  6a.  lustruiaented  Missile  Fiu  Test  Article:  View  A 


Figure  6b,  Instrumented  Missile  Fin  Test  Article:  View'  B 


Reproduced  From 
Best  Available  Copy 
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Figure  7a.  Missile  Fin  Heat  Pipe  Test  Setup:  View  A 


Figure  7b.  Missile  Fin  Heat  Pipe  Test  Setup:  View  B 
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Reproduced  From 
Best  Available  Copy 


Missile  Fin  Heat  Pipe 


1  DATE 

1  Oct.  1.96 

Oct  1.96 

Oct.  3,96 

Oct.  3.96 

Oct.8,96 

Oct8,96 

1  Oct.10.96 

Location 

■OOO 

579 

547 

546 

559 

520 

636 

345 

485 

2 

304 

344 

380 

370 

307 

295 

504 

154 

175 

3 

Cond.  RHP 

581 

608 

558 

556 

556 

552 

.  549 

213 

226 

4 

Evap  tube  top 

583 

584 

504 

503 

588 

581 

388 

537 

5 

Evap  tube  side 

636 

607 

517 

615 

650 

652 

#2HP 

112 

570 

6 

Evap  tube  bot 

596 

596 

578 

574 

586 

584 

551 

653 

7 

628 

602 

612 

610 

658 

661 

503 

624 

8 

Mid  tube  top 

470 

513 

516 

510 

501 

495 

242 

274 

9 

Mid  tube  side 

473 

510 

506 

495 

469 

460 

211 

260 

10 

Mid  tube  bot 

573 

584 

555 

549 

546 

541 

288 

294 

11 

Mid  tube  side 

518 

547 

520 

516 

503 

497 

273 

200 

12 

Cond  tube  top 

328 

366 

383 

375 

311 

302 

184 

160 

13 

Cond  tube  side 

392 

419 

428 

411 

353 

346 

166 

152 

14 

Cond  tube  bot 

500 

523 

495 

491 

485 

479 

152 

162 

15 

Cond  tube  side 

466 

494 

470 

466 

432 

426 

172 

160 

16 

574 

569 

573 

572 

591 

593 

332 

400 

17 

CondLE 

289 

331 

383 

371 

322 

311 

166 

156 

18 

636 

698 

639 

636 

701 

711 

723 

727 

19 

578 

570 

567 

565 

590 

590 

296 

414 

20 

Water  In 

25 

26.2 

24J2 

24.5 

24.4 

23.4 

23.4 

22 

21 

Water  out 

58.5 

59.5 

36 

35.6 

52.3 

49.4 

35 

38.7 

Flow  cc/min 

254 

254 

800 

908 

500 

563 

254 

254 

POWER 

596 

592 

661 

706 

976 

1025 

206 

297 

to  increase  the  likelihood  of  a  design  that  would  be  resistant  to  thermal  fatigue.  The  control  hole  tube/ 
envelope  assembly  was  formed  by  brazing  together  two  stainless  steel  pieces  with  wall  thickness  as 
small  as  0.012";  joining  was  accomplished  by  brazing  with  Nicrobraze  30.  Embrittlement  of  grain 
boundaries  in  the  braze  region  certainly  was  a  factor  in  the  wall  failures.  Thermacoie  constructed  heat 
pipes  using  wall  thickness  of  0.012"  or  thinner  twice  in  the  past.  One  design  used  a  laminated  hot 
isostatically  pressed  design  composed  of  alternating  layers  of  titanium  foil  and  a  mixture  of  titanium 
and  titanium  diboride  powder  (Rosenfeld  et  al.,  1995)  in  a  potassium  heat  pipe;  the  other  design  used 
a  tungsten  tube  in  a  lithium  heat  pipe  (Rosenfeld,  1992).  Both  of  these  approaches  have  been 
demonstrated  successfully.  It  is  believed  that  forming  the  missile  fin  envelope/control  hole  assembly 
from  a  single  electric  discharge  machined  (EDMed)  plate  would  be  a  better  fabrication  approach 
because  this  would  provide  a  single  wall  0.024"  thick  with  no  embrittling  braze  material  present  in  the 
wall.  Welding  rather  than  brazing  is  recommended  for  joining  the  end  plates  to  the  cover  plates.  Also, 
the  use  of  oxygen  getters  is  recommended  to  assist  in  suppressing  wall  grain  boimdaiy  failures. 

The  heat  pipes,  quartz  lamps,  power  supplies,  and  calorimeters  were  delivered  to  Wright 
Laboratory  for  further  testing.  This  test  report  is  being  supplied  to  document  the  test  results  and  design 
work,  to  assist  the  Air  Force  in  defining  future  work  in  this  area. 


6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

Two  stainless  steel  /  sodium  heat  pipe  missile  fins  were  designed,  built,  and  tested.  Each  heat 
pipe  operated  successfully  during  processing,  but  each  fin  &iled  in  a  similar  manner  during  subsequent 
testing.  Cause  of  failure  was  observed  to  be  a  failure  of  the  front  heat  pipe  envelope  wall  in  the 
location  near  the  control  hole,  vvdiere  the  wall  thickness  is  only  0.024". 

It  is  believed  that  by  modifying  several  design  parameters  and  fabrication  steps,  it  is  feasible 
to  construct  reliably  operating  heat  pipe  missile  fins  for  this  application.  A  summary  of  these 
recommended  changes  follows.  Thermacore  recommends  fabrication  of  additional  missile  fins  using 
the  following  changes. 
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1 .  Replace  the  stainless  steel  envelope  material  with  Haynes  230  superalloy.  This  material  has 
been  demonstrated  in  several  sodium  heat  applications  at  Theimacore  under  the  solar  dish  Stirling  heat 
pipe  development  program.  It  demonstrates  superior  weldability,  oxidation  resistance,  fatigue 
resistance,  grain  growth  resistance,  and  sodium  compatibility. 

2.  Constmct  the  missile  fin  assembly  from  a  single  block  of  metal  by  EDM,  rather  than  by  brazing 
together  thin-walled  subassemblies.  This  approach  will  provide  a  monolithic  wall  with  small,  crack- 
resistant  metal  grains.  This  approach  will  also  allow  construction  of  a  completely  welded  missile  fin 
assembly;  brazing  increases  the  likelihood  of  wall  embrittlement  for  thin-walled  assemblies  in  the 
braze  region. 

3 .  Use  a  piece  of  oxygen  getter  material  such  a  hafiiium,  zirconium,  or  titanium  in  each  heat  pipe. 
This  approach  will  reduce  the  likelihood  of  oxygen  attacking  grain  boundaries  and  crevices  inside  the 
heat  pipes. 
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Appendix  B: 
Fabrication  Drawings 
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Appedix  C: 

Methematica  Program  for  Transient  Model 


Front  Propag  (front  HP) 
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(*  This  workbook  calculates  the  transient  response 
of  a  Na  missile-fin  heat  pipe  based  a  lumped 
capacitance  model  of  continuTun  front  propagation. 


Front  Heat  Pipe  (leading  edge) . 


( *  constants  used  in  claculation  * ) 

tempvalue  s  {{O.O,  300.0}}; 
zvalue  s  {{0.0/  0.0}}; 

ze  -  0.050;  (*  evaporator  length,  m  *) 

zc  B  0.056;  (*  condenser  length,  m  *) 

thp  s  300.0;  (*  current  heat  pipe  ten^erature,  K  *) 

tmelt  =  371.0;  (*  melting  ten^erature  of  Na,  K  *) 

time  s  0.0;  (*  current  time,  s  *) 

dtime  s  0.1;  (*  time  increment,  s  *) 

ap  =  0.0001145;  (*  perimeter  by  thickness  of  pipe  wall,  m2  *} 
aw  =  0.000050;  (*  perimeter  by  thickness  of  wick,  m2  *) 

aes  0.0025;  (*  perimeter  by  evaporator  length,  m2  *) 

eps  -  0.55;  (*  porosity  of  the  wick,  non-dimensional  *) 

r  =  361.478;  (*  R,  gas  constant,  for  Na  *) 

gamma  —  1.67 t  (*  specific  heat  ratio  for  Na  *) 

1  B  0.015;  (*  characteristic  length  for  meem-free-path,  m  *) 

av  =  0.00020;  (*  area  of  vapor  section,  m2  *) 

he  =  752.0;  {*  convection  constant  for  evaporator,  W/m2-K  *) 

he  B  550.0;  (*  convection  constant  for  condenser,  W/m2-K  *} 

ter  B  1213.0;  {*  heat  source  ten^erature,  K  *) 

tinf  B  327.0;  {*  heat  sink  temperature,  K  *) 

p  B  0.05;  (*  perimeter  of  heat  pipe,  m2  *} 

zubO.O;  (*  initial  continuum  front  position,  m  *) 

hsl  B  113000.0;  (*  heat  of  fusion  for  Na,  J/kg  *) 


(*  temperature  dependant  functions  *} 
ge[t_]  :b  he  ae  (ter  -  t)  (*  W  *) 

Mg[t_]  :=  (5.335  (10^6)  -  4185.8  t  +  9.097  t^2  -  1. 1117* (10^-2)  t*3 
+  6.2049  (10*-6)  t*4  -  1.2942  (10*-9)  t*5)  (*  J/kg  *) 

rhoss[t_]  ;b  8055.0  (*  kg/m3  *) 

rhoni[t_]  :b  8900.0  (*  kg/m3  *) 

rhos[t_]  ;b  1033.45  -  0.221268  t  (*kg/m3  *) 

rholtt_]  :b  (991.42  -  8.8062  (10*-2)  t  -  3.4209  (10*-4)  t*2 
+  3.7279  (10*-7)  t^3  -  1.9426*(10*-10)  t*4 
+  3.8804  (10*-19)  t*5)  (♦  kg/m3  *) 

rhov[t_]  :b  E3cp[-86.67  +  0.3067  t  -  4 . 972*  (10*-4)  t*2 

+  4.323*(10*-7)  t*3  -  1.925*(10*-10)  t*4 
+  3.443*(10*-14)  t*5]  (*  kg/m3  *) 

cpss[t_]  ;b  372.115  +  0.421348  t  -  0.000188779  t*2  (*J/kg-K  *) 

cpni[t_]  :b  235.0  +  0.863643  t  -  0.000551786  t*2  (*J/kg-K  *) 
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cpsodl[t_]  :=  1550.0  -  0.564899  t  +  0.000274368  t*2  {*J/kg-K  *) 

cpsods[t_]  :=  1364.2  (*  J/kg-K  *) 

;=  (6.8708  (lO^^-O)  +  2.506  (lO^^-O)  t  -  3.6443  (10^-11)  t^2 
+  4.9847  (10^-14)  t^3  -  3.217  (10^-17)  t^4 
+  7.6982  {10'‘-21)  t^5)  (*  N-s/ib2  *) 

qs[t_]  :=  av  rhov[t]  hfg[t]  Sgrt[0.5  gamma  r  t/ (gamma  +  1)]  (*  W  *) 

cl[t_]  :=  (rhoss[thp]  cpss[thp]  ap  +  (1.  -  eps)  (rhoni[thp]  cpni[thp]  aw) 
+  eps  rhos[th.p]  cpsods[thp]  aw)  (*  J/m-K  *) 

c2 [t_]  :=  (rhoss[thp]  cpss[thp]  ap  +  (1.  -  eps)  (rhoni[thp]  cpni[thp]  aw) 
+  eps  rhol[thp]  cpsodl[thp]  aw)  (*  J/m-K  *) 

(*  stage  1  *) 

While [  (*  test  *)  thp  <  tmelt, 

(*  body  *) 

thpnew  =  thp  +  (ge[thp]  dtime) / (cl [thp]  ze) ; 

thp  s  thpnew; 

time  s  time  +  dtime; 

AppendTo[tempvalue,  {time,  thp}]; 

AppendTo[zvalue,  {time,  zn}]; 

(*  end  of  While  loop  *)  ]; 

Print [time,  ”  s  for  the  evaporator  to  reach  melt  temp."] 

Print [thp,  "  K  is  the  current  evaporator  tesqperature . " ] 

(*  end  of  stage  1  *) 

1.3  s  for  the  evaporator  to  reach  melt  temp. 

373.76  K  is  the  current  evaporator  temperature. 

(*  stage  2  -  melting  evaporator  *) 
dtimemelt  -  eps  rhol [500.0]  aw  ze  hsl/ge[thp]; 
time  =  time  +  dtimemelt; 

Print [time,  "  s  to  thaw  evaporator  section."] 

Print [thp,  "  K  is  the  current  evaporator  teng>erature . " ] 

AppendTo[teiig>value,  {time,  thp}]; 

AppendTo [zvalue,  {time,  zn}]; 

(*  end  of  stage  2  *) 

1.38827  s  to  thaw  evaporator  section. 

373.76  K  is  the  current  evaporator  temperature. 
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(*  stage  3  *) 

(*  first  find  transition  temperature,  tt  *) 

tt  =  X  /.  FindRoot[x  ==  5020.0  Pi  (muv[x] / (rhov[x]  1))^2.0, 

{X,  600.0}]; 


(*  iterate  untill  tt  is  reached  *) 

While [  (*  test  *)  thp  <  tt, 

{*  body  *) 

thpnew  =  thp  +  (ge[thp]  dtime) / (c2 [thp]  ze); 

thp  s  thpnew; 

time  =  time  +  dtime; 

AppendTo[teinpvalue,  {time,  thp}]; 

AppendTo[z value,  {time,  zn}]; 

(*  end  of  While  loop  *)  ]; 

Print [tt,  "  K  is  the  transition  temperature."] 

Print [time,  "  s  to  reach  transition  temperature."] 

Print [thp,  "  K  is  the  current  evaporator  temperature."] 

(*  end  of  stage  3  *) 

832.091  K  is  the  transition  temperature. 

15.1883  s  to  reach  transition  temperature. 

832.717  K  is  the  current  evaporator  temperature. 
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{*  stage  4  *) 

dtime  =  10.0  dtime; 

zn  s  ze; 
znew  =  ze; 
zold  s  0.0; 

(*  iterate  untill  zn  s  ze+zc  or  steady-state  is  reached  *) 


While[  (*  test  *)  znew  <  (ze  zc)  &&  Abs[zn  -  zold]  >  0.000001, 

{*  body  *) 

{thpnew,  qf}  =  {xx,  yy}  /.  FindKoot[ 

{qe[thp]  ==  he  p  (zn  -  ze)  (thp  -  tinf)  +  yy 
+  c2 [thp]  zn  (xx  -  thp) /dtime, 
yy  ==  Min[qe[thp] ,  qs[thp]]  -  he  p  (zn  -  ze)  (thp  -  tinf) 
-  c2[thp]  (zn  -  ze)  (xx  -  thp) /dtime}, 

{xx,  800.0},  {yy,  999.0}];  (*  K,  and  W  *) 

znew  =  zn  +  qf  dtime/ (rhol [500.0]  eps  aw  hsl  +  cl [350.0]  (tmelt 
-  tinf)  +  c2 [700.0]  (thpnew  -  tmelt));  (*  m  *) 

thp  =  thpnew; 

zold  B  zn; 

zn  =  znew; 

If[Abs[zn  -  zold]  <=  0.000001,  Print ["steady-state  reached  at"]]; 

time  s  time  +  dtime; 

AppendTo[tei^pvalue,  {time,  thp}]; 

AppendTo[z value,  {time,  zn}]; 

(*  end  of  While  loop  *)  ]; 


Print [time,  "  s  after  startup."] 

Print  [thp,  "  K  is  the  heat  pipe  temperature."] 
Print [zn,  "  m  is  the  continuum  front  location."] 

(*  end  stage  4  *) 


steady-state  reached  at 
205.188  s  after  startup. 

832.717  K  is  the  heat  pipe  temperature. 
0.101383  m  is  the  continu\am  font  location. 
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(*  stage  5  *) 

(*  iterate  untill  steady  state  is  reached  *) 

thpnew  -  tinf; 
thpold  =  tinf  -  10.0; 

While [  (*  test  *)  N[thpold,  4]  !=  N[thpnew,  4], 

(*  body  *) 

thpnew  =  thp  +  (Qe[thp]  -  he  p  zc  (thp  -  tinf))  dtime/ (c2[thp]  (ze  +  zc))i 

thpold  =  thp; 

thp  =  thpnew; 

time  =  time  +  dtime; 

AppendTo[teinpvalue,  {time,  thp}]; 

AppendTotz value,  {time,  zn}]; 

(*  end  of  While  loop  *)  ]; 

Print ["Steady-state  reached  at"] 

Print [time,  "  s  after  startup."] 

Print [thp,  "  K  is  the  heat  pipe  temperature."] 

(*  end  stage  5  *) 


1. 

Steady-state  reached  at 
585.488  s  after  startup. 

814.041  K  is  the  heat  pipe  temperature. 


ListPlot [Take [ten^value, 250] ] 
ListPlot [Take [zvalue,  250] ] 
IiistPlot  [tempvalue] 

ListPlot [zvalue] 
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